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Abstract 

In e + e~ collisions the particle spectra at low momenta reflect the properties of the 
underlying "soft" QCD gluon bremsstrahlung: the particle density, in the limit p — > 0, 
becomes independent of the incoming energy yfs and directly proportional to the colour 
factors Ca, Cp for primary gluons or quarks respectively. We find that experimental data 
from the pp and nuclear reactions reveal the same behaviour: in the limit pt — > the 
invariant particle spectra become independent of the collision energy, and their intensities 
in e + e~, pp and nuclear reactions are compatible with the expected colour factors Cp: 
Ca'- (Npart/tyCA f° r N par t nucleons, participating in the interaction. Coherent soft gluon 
bremsstrahlung is, therefore, suggested to be the dominant QCD mechanism for the soft 
particle production in all these reactions. These "soft" particles probe the very early 
stage of hadron formation in the collision. Future measurements at the LHC will provide 
crucial tests on the contributions from possible incoherent multi-component processes. 

1 Introduction 

One of the interesting phenomena in high energy collisions is the production of soft particles 
with low momenta. It is usually described in the existing Monte Carlo models by a non- 
perturbative mechanism which transforms the partonic final state into a corresponding final 
state of hadrons. Best understood is the e + e~ process where the primarily produced qq pair 



evolves into a partonic cascade and finally into hadrons. A similar description applies to the 
hadronic collisions with the high p? jets, where the partonic cascades are initiated by the 
active partons and the parton spectators; also an incoherent superposition of such processes 
could become important at high energies in the TeV region. In nuclear collisions the collective 
phenomena play an important role as well. 

This paper is focused on the universal feature of the soft (p —> 0) particle production in 
all these collision processes which is based on the universality of the soft gluon bremsstrahlung 
with the intensity depending on the directions and colour charges of the participating partons 
but not on the collision energy. From this point of view the perturbative analysis of the quark 
gluon cascade (see, for example, [TJ |2]) is naturally combined with the idea of a close similarity 
between the parton and hadron spectra ("Local Parton Hadron Duality" - LPHD |3J). Such 
a scenario is deeply rooted in the space-time picture of the parton-cascade development and 
hadronization, see [TJ 0]. This approach has been developed and tested in details for the 
hadronic final states in e + e~ annihilation verifying, in particular, the energy independence of 
the soft limit and the expected dependence on the primary colour charges (quark or gluon jets), 
but it should also be applicable to other hard processes including ep and pp (pp) collisions with 
high p T particle production, where the soft particles belong to the "underlying event" [5] - [TJ. 
A distinct dependence on the primary colour charges makes this approach quite different from 
the simple process-independent universality, for example, of a thermodynamical origin. 

Here, we are extending this approach to the description of the so-called "minimum bias" 
events in the hadronic and nuclear collisions where we expect similar effects generated by the 
primary soft gluon bremsstrahlung. Note, first, that it is plausible to expect that at high 
energies the bulk of the soft particles originates in a semihard process from the fragmentation 
of the so-called minijets, which are produced by the perturbative mechanism similar to the 
e + e~ case. Next, recall that the perturbative dynamics results in the limiting behaviour of 
the soft particle densities, since the radiation of soft gluons from several quarks and gluons 
develops in a coherent way: different sources of secondaries should act as a single source with 
an effective colour 'charge' equal to the vector sum of all the colour 'charges', since a gluon 
with a large wave-length cannot resolve the smaller details (individual elements of the 'colour 
antenna'). It is worth emphasizing also that the low momentum hadrons are formed at the 
initial stage, before the formation of fast particles [TJ 0] . Thus, the low-p hadrons emitted at 
small rapidities in the centre-of-mass frame (c.m.s.) bear no information about the complicated 
structure of the whole event, in particular, about the secondaries in the fragmentation regions 
and the properties of the initial hadrons. 

It is natural to expect the analogous interference phenomena in the light hadron (pion) 
emission, when the particle inverse momentum, 1/p, which is a characteristics of the resolution, 
exceeds the size of the interaction region. Thus, we should observe the limiting behaviour in 
the soft particle densities in the proton-proton, proton-nuclear and nuclear-nuclear collisions 
similar to that observed in the e + e~ annihilation. All these processes can be characterized by 
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the limiting density I of the soft particle spectra [5] at vanishing rapidity y and pt 



norm 
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In the case of e + e~ annihilation the normalization cross section is cr toi , while for the pp and 
nuclear collisions the rates are often normalized by Oi n = a to t — &eh However, for the studies of 
soft particle production it is more appropriate to use for normalization another quantity, namely 
the so-called non-diffractive cross section, a^D, defined as (Tnd= °in — &d- Here <7d describes 
the corresponding contribution from the single diffractive and double diffractive dissociation. 
Strictly speaking, on depends on the collision energy (see, for example, a compilation in Ref.[8] 
and the analysis in [9]). However, this energy variation is within the uncertainties of the current 
data, and it is plausible to evaluate the non-diffractive cross section <Jnd as being 10% lower 
than £Tj n . Note also, that here and in what follows the quantity 1$ is defined as the density 
averaged over the charged particles [\{h + + h~)]. After recalling the relevant observations for 
the e + e~ collisions we will turn to the study of the limiting behaviour of the density ([1]) in the 
pp and AA collisions. We assume that these processes are initiated by the gluon exchanges, 
but more general configurations are possible. 

The relation between parton and hadron final states is formulated at first for inclusive 
observables such as we discuss here. It is an interesting question how such relations can be 
realised by the exclusive final states including resonances. In the limit p — > soft gluons in the 
cascade add coherently and are represented by the Born term; likewise the soft pions from a 
general final state should add coherently to yield the energy independent density. 

Analogies between the bremsstrahlung in QED and multiple particle production in hadronic 
collisions date back to the pre-QCD times [TUl [11] and has lead in particular to the concept 
of flat rapidity plateau. Within QCD such an analogy was further developed [12] and the 
picture with one-gluon exchange initiating high energy pp collisions has appeared [13]. The 
present QCD-based approach relates hadrons to the full quark-gluon cascade which generates 
particle jets, and the "classical" bremsstrahlung Born term is dominant only in the ultrasoft 
limit pt — > 0. 

2 Lessons from the soft particle spectra in e + e~ collisions 

Let us start from recalling the behaviour of the momentum spectrum of particles in e + e~ 
annihilation in the soft limit as the benchmark reaction. In this process it is convenient to 
measure the inclusive momentum spectrum since this does not require a determination of the 
jet axis which is needed for the spectra defined in terms of y and pt- 

1 We refer here to the one-side limit y — > 0, not to the double-side limit \y\ —> 
2 In some cases, when we refer to the charged pions, the value of Iq is 7% lower. 
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The existing data can be well described within the Modified Leading Logarithmic Approxi- 
mation (MLLA) for the QCD parton cascade together with LPHD for the description of hadrons 
(see, for instance [HE])- For the very low momenta a model dependent prescription is applied 
allowing to account for the mass effects. The perturbative expansion of the momentum spec- 
trum of partons can be derived by an iteration of the Born term in the respective evolution 
equation. The first term in MLLA has the form [5], H] 



1 - C A,F , 2 ( 2 ) 



while the higher a s corrections, which depend on the primary energy, vanish as ln(ln(fcr/A)/ 
ln(<5o/A)) when the transverse momentum of the gluon, kx, aproaches its minimum (cut-off) 
value Qq. Here n^' F is the multiplicity of gluons emitted by the configuration with the gluon (A) 
or quark (F) colour charge; y = ln(l/x), momentum fraction carried by the gluon is denoted as 
x and QCD scale - as A. In this limit the Born term remains and it is independent of primary 
energy and directly proportional to the colour factors. These features stay unchanged if the 
spectra are slightly modified at low momentum within LPHD to include particles with mass 
and one arrives at the perturbative expectation for the limiting behaviour [5j [6] 

dn , . 

p -> : E— > const. (3) 

d 6 p 

In this limit the gluon has a large wave-length and only resolves the primarily produced partons, 
a radiation process represented by the Born term. Another feature of the approximation (j2j) is 
the "flat rapidity plateau" , that is, the independence at fixed k T of the particle density on y or 
angular rapidity y. This approximation is true as far as we neglect the kinematical (threshold) 
effects, when at large y the energy of the observed secondary particle becomes comparable with 
the energy of the initial interaction. 

The available data show, that at the lowest accessible particle momenta, around p = 
0.2 GeV, the invariant density E-^ for pions rises only by about 20% over the c.m.s. en- 
ergy range ^/s = 3.0 -r- 160 GeV, and with increasing momentum p this rise becomes stronger, 
see [5] - [TJ. The MLLA formulae provide a rather good fit of the total energy and momen- 
tum dependence of the observed one-particle spectrum, and clearly demonstrate the approach 
towards energy independence of the spectrum as in the limit (J3j). Therefore, the observed rise 
of the particle density at mid-rapidities and the rise of the overall multiplicity with energy 
increasing is caused by the high px particles. Note that in the soft limit p — y one obtains 
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y-*0,PT— >0 



p-s-0 

since the invariant particle distribution approaches a constant value. 

The second important lesson from the e + e~ data concerns the dependence of the soft particle 
density on the primary colour charges, that is, on the difference between the primary quark- 
antiquark and the gluon-gluon colour antennae-dipoles. This was confirmed by the study of 
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3-jet (qqg) events at different inter-jet angles. In the limit when the gluon is parallel to a quark 
(antiquark) the soft radiation pattern is the same as that in the case of a qq dipole, since the 
soft gluon cannot resolve the parallel quark and gluon. On the other hand, if the qq system 
recoils against the gluon, it is seen by the emitted soft gluon as a colour octet source. Then 
the soft radiation density in these two limiting cases should differ by the ratio of colour factors 

with Ca/Cf = 9/4, in agreement with (j2]). Experimentally, one cannot reach these extreme 
limits but we can derive predictions on how the soft radiation density perpendicular to the event 
plane in the 3-jet event varies with the inter-jet angles between the two extremes [5]. These 
expectations are well confirmed by the DELPHI data [14J which successfully reproduce the ratio 
(JSJ). Note that the corresponding ratio for the global quark and gluon jet multiplicities strongly 
deviates from ([5]) because of the influence of the particles with higher and the importance 
of higher orders in the perturbative MLLA calculation. 

Neglecting the effects of order of 1/Nq, the primary qqg antenna can be represented by the 
two qq dipoles, where the gluon is replaced by the parallel q and q. In this approximation the 
ratio ([5]) becomes Ca/Cf — > 2, corresponding to the number of radiating dipoles. The radiation 
of soft gluons perpendicular to the dipole varies with the opening angle as Ejj§^ ~ (1 — cos 0). 
The experimental results demonstrate that the soft particle density follows the lowest order 
QCD expectations for the given configuration of colour charges, see for a review [I]. 



3 Soft particle spectra in high energy pp (pp) collisions 
3.1 Minimal model for the minimum bias events 

In analogy to the e + e~ process, let us start from the minimal partonic process, which can be 
responsible for the very soft gluon bremsstrahlung. We assume here that the underlying physics 
of the minimum bias events is based on the collisions of two partons within the protons. The 
partonic process of lowest perturbative order corresponds to one gluon exchange, which leads to 
a dominantly small scattering angle and a non-vanishing cross section at high energies. In the 
case of elastic scattering between the two incoming quarks the exchange of the t-channel gluon 
rearranges the incoming colours and creates the outgoing colour charges which leads to the 
radiation of the soft gluons from the effective colour octet dipole. For a large number of colours 
in the above-mentioned approximation this radiation can be described as being generated by a 
superposition of the two aligned dipoles. In the case of small-angle scattering the same soft 
gluon radiation pattern appears also with the incoming quarks replaced by gluons, as can be 
explicitly seen from the radiation patterns given in [15] . 

Therefore, the limiting soft radiation densities in the pp and e + e~ collisions should differ by 
a factor of 

p : I% p /lt e ~ ~ Ca/Cf (6) 
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at all energies, as in the ratio in Eq. (E]) for the spectra induced by the gg and qq dipoles. 

For a similar reason the particle multiplicity ratio for the two processes is expected to be 
equal to the same quantity Ca/Cf [IS]- In our case, we anticipate that this QCD expectation 
could be valid only in a specific soft limiting case, while as we discuss, the pr-behaviour of the 
spectra in both processes are different and energy dependent, therefore the total multiplicities 
should differ as well. 

One can also consider more complex situations corresponding to the multiple gluon ex- 
change. However, if there is a multiple gluon exchange between the pair of quarks, then the 
total colour of the exchanged system could be only an octet or a singlet; since the singlet ex- 
change rather corresponds to a diffractive process, its contribution to the central production 
is negligible, and we come back to the previous case. If the multiple gluon exchange involves 
both quarks and spectator diquarks, but does not destroy the diquark, (which acts as a local 
object), then the colour of the 'quark-diquark' system can be either an octet or a singlet, and, 
thus, again after the multiple gluon exchange this generates in the t-channel only the 'octet' 
colour flow. Finally, recall that in the case of more complicated diagrams the dominant Leading 
Logarithmic contribution in both the DGLAP and BFKL evolutions also corresponds to the 
colour octet exchange. 

In the early discussions [T7] several specific models were proposed, which predicted the ratio 
of the central rapidity particle densities in the pp and e + e~ collisions. In these models the final 
state in pp-interactions is constructed from the two chains of particles connecting q and qq 
while there is one chain between q and q in e + e~ collisions, but the ratio of these densities 
in both processes depends on the amount of overlap between both of the chains and can vary 
between one and two. Later on, the semi-hard gluon bremsstrahlung has been included in the 
models and this leads to a rapidity plateau rising with energy, and a simple relation between 
both processes fades away. An application of the Lund model to hadronic scattering at the 
energies below ISR was discussed in [18], where the same one-string description was used for 
both, e + e~ and pp collisions. 

At higher energies this scenario should be modified in order to incorporate the multiple 
interactions (MI) of partons (that is simultaneous interactions of two or more pairs of par- 
tons) and minijet formation, which would lead to a certain energy growth of Io f|, for recent 
development based on the dipole cascades and multiple collisions, see 



Current Monte Carlo models include the MI option, which in terms of a simple eikonal 
model for soft pp-scattering, where one interaction corresponds to the exchange of a cut bare 
Pomeron, is described by the contribution of the few cut Pomerons. Particle density produced 
in the central region by one cut Pomeron is energy independent @. Recall, that the cross section 



3 We thank Gosta Gustafson for a discussion of the Lund model results. 

4 In a simple soft scattering model particle density is energy-independent. For shower MCs, where each 
MI is modeled by the perturbative QCD contribution of the form of / dx\dx2g{x\, qt)g(x2, Qtjdfr/d qt \q t >q min , 
the approximately constant behaviour results from a compensation between the growth of the gluon density 
xg(x,q t ) ~ x~ x with decreasing momentum fraction x ~ qt/y/s and an increase of a cutoff q m i n ~ l/s°- 08 ~ 013 . 
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of hard subprocess is da/d 2 q t oc 1/qf- However, in the Gribov-Regge theory the mean number 
of Mis, that is the number of cut Pomerons, increases as s A , where the bare Pomeron intercept 
is «p(0) = 1 + A (in terms of the gluon density xg(x, q t ) ~ x~ x we identify A = A) @. This is 
the well known AGK result [21] - single particle inclusive cross section is described just by one 
Pomeron exchange, and, thanks to the AGK cancellation, there is no absorptive corrections. 
Therefore, we arrive at da/d 3 p oc s A , while the growth of inelastic cross section is reduced by 
the absorptive effects (asymptotically instead of Oi n oc s A we should reach the Froisart limit 
(Ji n oc In 2 s). Phenomenologically, the energy behaviour of Oi n is close to s t with e ~ 0.08 -j- 0.1. 
Thus, the particle density Iq = {1/ ai n )Eda / d?p oc s A ~ e , and this growth with energy reflects 
the increasing number of Mis, which act in the present Monte Carlo models incoherently. An 
account of the coherence effect should diminish the multiplicity of low pr-secondaries produced 
by MI. An explicit computation [22] using Perugia P0 tune [23] of the Pythia MC model, shows 
that in the c.m.s. energy interval y/s = 23 -r- 14000 GeV the invariant density in the pp collisions 
E-^ in the range pt = 4- 0.5 GeV and y = -4 0.5 is rising by a factor of about 2. Other, 
'more aggressive', tunes may result in even higher rise, up to a factor of 3 at most B As we 
discussed above, without incorporation of coherence, such increase should hold on up to the 
asymptotic energies. 

Finally, from such simple properties as eq. (EJ) one should not conclude that these two 
processes are intrinsically very similar beyond the considered limit pt — > and, indeed, the 
above models show very different phenomena. In fact, the event structure of a 2-jet system in 
e + e~ annihilation is not the same as that of a system with two (q-qq) pairs in pp collisions with 
its much larger fluctuation properties. However, in the considered limit, the produced partons 
(q-qq) are resolved as an extended colour octet source with the enhanced radiation according 
to eq. (jUJ). Such intuitive limits may help to understand the low p T production based on a 
perturbative ansatz. 



3.2 Energy dependence of the limiting soft spectrum 

For pp interactions in case of the minimal model discussed above the soft radiation pattern 
should again be energy independent analogously to the e + e~ case, but with the intensity about 
a factor of two higher. We are analyzing below the p^-spectra in pp collisions in order to 
determine the quantity Iq for px — >• at the central rapidity y — 0, as defined by relation 
([1]). The previous analysis of the invariant cross section m the p*r-range 0.3 4 1.0 GeV 

indicated a moderate rise with energy increasing from the ISR to the SppS colliders, see Ref. 
[7]. Later on, further data from CDF at the Tevatron at higher energies up to 1800 GeV 
and from RHIC became available, which required a fresh look at this problem. 



5 In the recent model for soft hadronic interactions [9] A ~ 0.3, which is in agreement with the resummed 
NLL BFKL result, for example [20] . 

6 We are very grateful to Torbjorn Sjostrand for the discussion of this and other issues related to results in 
Pythia model. 
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The inclusive spectra were measured down to the transverse momenta of 0.1 -j- 0.4 GeV. 
Some experimental groups performed fits using certain simple parametrizations of the spectra 
from which the quantity Iq can be easily determined. The British Scandinavian Collaboration 
(BS) [25J fitted the invariant cross sections at different energies between y/s = 23 to 63 GeV as 

E^- = Aexp (Bp T + Cp% + Dy 2 ), (7) 



while the UAl [26] and CDF [24] groups used the parametrization 



E^ = A(l+p T / Po r n . (8) 

From these fits the quantity ([I]) can be found as Iq = A/<r in . The BS collaboration |25j 
confirmed the exponential behaviour down to the pion momenta as low as pt ~ 0.1 GeV. 
Correspondingly we derived the quantity J from the STAR data |27j using the exponential 
extrapolation of their charged particle spectra. 

As the results depend on the function chosen for extrapolation we also considered other 
functional forms. The PHOBOS [39] collaboration has measured the spectra down to very 
small pt = 0.03 GeV in AuAu collisions at 200 GeV. These data indicate a flattening of the 
distribution towards small Pt according to a functional form following from a thermal model 

„ dn A 

E = (9) 

d 3 p exp(mr/T) — 1 



with ttit = \j m 2 + p\ and temperature T as parameter. 

In order to test the dependence of the extrapolated results on the choice of the function we 
compared the pp data of highest accuracy which is provided by the STAR collaboration [27] 
with errors of ±3% to an exponential as in ((7j) (but with C = 0) and the function fl9]) in the 
Pt ranges 0.225-0.425 GeV and 0.225-0.625 GeV. We found that the extrapolated cross section 
at pt = is reduced by 25% for the thermal fit. 

It should be noted that the forms ([TJ and (jSJ) do not have the correct analytic behaviour 
at pt — > where they behave linearly in pt and the invariant cross section is not regular; 
therefore, one expects to obtain an overestimate from such fits. This problem does not occur 
with the non-singular fits of the type ([9]) with quadratic behaviour in pt- 

Experimental results over a range of energies are only available from the exponential and 
power like fits which we consider next. To obtain the cross section, <7j n , we take the total and 
elastic cross section data for the pp (BS) and pp scattering (UAl, CDF), as collected by the 
Particle Data Group [28]. The results are shown in Table 1 after an interpolation. In order to 
determine the quantity Iq for the STAR entry in this Table, we renormalized their minimum 
bias data to the a in as well. 

The data for Iq in the last column of Table 1, shown in Fig. (TJ indicate an initial rise in the 
BS data with increasing energy, but it does not hold at higher energies. Remarkably, the data 
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Table 1: Soft limit Iq = Ajoi n from the exponential fits to the single charged-particle [(h + + 
h~)/2] spectra in the pp (pp) collisions; for the BS and STAR data we sum over the 7r ± , K ± ,p ± 
distributions. 



Exp 




PT,min 


A 


Otot 


0~el 




Iq = Ajo in 




[GeV] 


[GeV] 


[mb/GeV 2 ] 


[mb] 


[mb] 


[mb] 


[GeV~ 2 )] 


BS 


23 


0.1 


191 ±7 


39.4 


6.8 


32.6 ±0.5 


5.9 ±0.3 


BS 


45 


0.3 


238 ±7 


41.9 


7.5 


34.4 ±0.7 


6.9 ±0.3 


BS 


63 


0.1 


307 ±20 


43.0 


7.8 


35.2 ±0.6 


8.7 ±0.7 


STAR 


200 


0.2 










7.5 ±0.8 


UA1 


200 


0.25 


286 ± 17 


52 


9.2 


43 ±4 


6.6 ±0.7 


UA1 


500 


0.25 


408 ± 24 


62 


13 


49 ±2 


8.3 ±0.6 


CDF 


630 


0.4 


300 ± 20 


63 


13 


50 ±2 


6.0 ±0.5 


UA1 


900 


0.25 


382 ± 20 


68 


15 


53 ±4 


7.2 ±0.7 


CDF 


1800 


0.4 


450 ± 10 


74 


17 


57 ±3 


7.9 ±0.5 



in the large energy range of 40 -j- 1800 GeV show a rather energy-independent behaviour, and 
they fluctuate around some common mean value by about ±15%: 

inelastic pp/pp collisions (exp. fit) : J ~ (7 ± 1) GeV~ 2 (10) 

and for minimum bias events with the above estimate for diffractive cross sections 

non — diffractive pp/pp collisions (exp. fit) : Iq ~ (8 ± 1) GeV~ 2 (11) 

Next we consider the extrapolation with the thermal function (Q. While our main finding, the 
lack of any major energy dependence of Iq does not depend on the extrapolation procedure the 
central value of Iq will be reduced by 25% and we obtain 

non — diffractive pp/pp collisions (therm, fit) : Iq ~ (6 ± 1) GeV~ 2 (12) 

This observation of energy independence is similar to that in e + e~ annihilation, discussed 
above, and it is consistent with the expectation corresponding to a single coherent bremsstrahlung 
process. On the other hand, for an incoherent superposition of many processes one would rather 
expect Iq values, which rise with energy increasing since the number of such processes would 
typically increase. We should also recall that the very weak energy dependence (if any) of the 
quantity Iq contrasts to the energy dependence of the pr-integrated central rapidity density 
dn/df]\ v= Q. The latter observable is rising in the considered energy range by a factor 2-3 in the 
pp and AA collisions (see, for example, the recent discussions in [29]). 

Since the soft gluon emission is driven by the lowest order diagram, according to eq. ([2]), 
we may expect that the same limiting behaviour (da / dyd 2 pT — > const for p T — > 0) should be 
observed at any rapidity and not only at y = and we expect a flat rapidity plateau in this 
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10 100 1000 

[GeV] 

Figure 1: Soft limit Iq of the invariant density i?Jr of charged particles [(h + + h~)/2] in pp 
collisions as a function of c.m.s. energy y/s (from exponential extrapolation, Table 1). 

limit, as mentioned above. This is true, as long as we can neglect the kinematical (threshold) 
effects for large rapidity y. Also, we are looking for the soft hadrons (pions) which are formed 
first and do not participate in further collisions with the faster secondaries. These conditions 
are better satisfied for the particles in the central region. For our analysis we use the inclusive 
cross sections at y = where we have the best collection of data. Actually, the ^-dependence 
was studied by the BS collaboration [22] at the ISR. In the rapidity range \y\ < 1.5 and pr < 1 
GeV this dependence is rather weak and the Gaussian fit as in eq. (7) yields values around 
D ~ 0.1. Yet smaller values of D are required if additional data from large y are included. 
Measurements at higher energies would be of interest as well. 

Finally, we comment on the first publication of charged particle pr spectra in non-diffractive 
minimum bias events from the LHC at \fs = 0.9 and 2.36 TeV by the CMS collaboration [30J. 
One can observe the convergence of the spectra for px — > 0. An exponential extrapolation of 
the data to this limit yields Iq ~ 9 GeV -2 , which is consistent with (ITT]) . Further studies at 
higher energies at LHC will provide a critical test of the contributions of additional incoherent 
multiple interaction components. 

3.3 Comparison with the e + e data 

It is also instructive to compare the observation ( fTTl) with the corresponding quantity for the 
e + e~ annihilation. As we already discussed, since in the pp collisions the diffractive events 
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Figure 2: The density 1(E) = j^p^ of pions [(tt + + n~)/2] as a function of the c.m.s. pion 
energy E, as derived from the fit to the BS pp data and normalized to the non-diffractive cross 
section at 53 GeV. 

lead only to a small contribution to the central rapidity region, the diffractive piece should be 
subtracted in this comparison. We estimate the limiting value I$ +e in three different ways: 

1. The TPC/27 collaboration [31] compared the results for inclusive n/K/p particle spectra, 
measured in the e + e~ annihilation at 29 GeV, with the corresponding spectra from the non- 
diffractive pp collisions at 53 GeV, obtained by the BS collaboration [25] . The higher pp collision 
energy was chosen in order to take into account the lower effective energies of the parton-parton 
collisions. 

The spectra measured in both reactions have comparable magnitude in the transverse mo- 
mentum range 0.25 < pr < 0.5 GeV, where the rapidity y and px in e + e~ annihilation events 
are defined with respect to the sphericity axis. As discussed above, the pp data [25] follow the 
exponential behaviour for small p? < 1 GeV down to pt ~ 0.1 GeV, and one can extrapolate 
the data for pions (7r + or n~ reproduced in [31] )@ to I ~ 8 GeV~ 2 in agreement with (TIT]) or 
to a 25% lower value if we take into account mass effects as above. On the other hand, the 
PT-spectra in e + e~ annihilation [31] deviate from the exponential behaviour, decreasing more 
weakly for higher p? and flattening towards the small p^-values with an extrapolated limit 
Iq ~ (3.0 ± 0.3) GeV -2 . This way we arrive for the pions at the ratio r = Iq P /Iq +£ 2.7 or at 
r rs 2.0 if we apply again the 25% reduction for mass effects in the pp analysis. For the kaons 
an exponential extrapolation appears to be applicable for both processes at px < 1 GeV, but 
with different slopes, and the corresponding values are I ~ 0.48 GeV -2 and I ~ 0.24 GeV~ 2 
so that Iq^ /Iq e ~ 2 in the case of kaons as well. While the pion and kaon spectra are steeper 
in pp collisions and cross those in e + e~ collisions, the two spectra for protons do not show this 
behaviour but are rather similar in the considered range. Then we look at other results as well. 

2. Since no other group presented their results in terms of the y,PT variables in the e + e~ 
experiments, we consider here also the extrapolation in the standard variables, momentum p or 

7 The data plotted in Fig. 2 of [3T] should be divided by 2tt in order to obtain the quantity —gfe- 
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the energy E, which are directly measurable without any reference to a jet axis. Then in the 
e + e~ collisions one considers the distribution t-^tt 1 . Such spectrum can be also determined 

inp z dp L 

in the pp collisions, and we calculate its shape using the parametrisation for the invariant 
spectrum of the pions [(/T + + n~)/2] at 53 GeV Ejjg^ = ^ d ^ p2 = I(y,pr) cr, as given by the 
BS collaboration [25], see Eq. ([I]), with the parameters A = 212 mb/GeV 2 , B = -7.3 GeV" 1 , 
C = 1.2 GeV -2 and D = —0.13 and with normalization by the non-diffractive cross section 
as in [21]. We calculate the corresponding distribution over the c.m.s. momentum p, using 
dy = dp z /E as 

f P = 2 I**I *^*-v^> < 13 > 

E dn f p dp 



Aixp 2 dp 



[ —i(y(p>Pz),PT(p,Pz)) (14) 

Jo P 

where y = ~ In and pt = \/p 2 — pi- This spectrum is shown in Fig. |2] as a function of 

c.m.s. pion energy E. One can see the peak structure at low energy E q Pion distributions 
over the c.m.s. energy, E, from the e + e~ annihilation are collected in Fig. 3a of [5]. The 
available data, especially ARGUS results [32] at y/s ~ 10 GeV, extending towards the small 
momenta p ~ 0.05 GeV, do not show this kind of peak. They rather show a flatter exponential 
distribution over E with the limit Iq ~ 3 GeV -2 for E — > m^, about half the pp value. § 
The e + e~ data at higher energies E are rising with */s and will cross eventually the spectrum 
for pp collisions in Fig. [2j Note that the dependence on the jet-axis definition disappears for 
p — > 0. These results demonstrate again that the inclusive spectra in e + e~- are flatter than in 
pp collisions at the low values of E. 

3. The TASSO collaboration [33J has presented fits to the invariant spectrum with expo- 
nential form in E 

JLf^A^M-BnE), (15) 

m 

with 2 or 3 terms, which we can use to obtain Iq from the limit E — > m n . For normalization 
we take the total cross section as a to t = with R 4 representing the data collection in 

[28] in this energy region. 

Results based on the TASSO fits and our estimates from the experiments at lower e + e~ 
energies with small momentum cut-off are summarized in Table [2j They show the consistent 
values for Iq. At higher energies no fits allowing to perform the extrapolation are available. 
However, as already noted in Sect. 2, the energy dependence at p — 0.2 GeV is weak over 
the energy range yfs ~ 3 4- 160 GeV, and the deviations from the model incorporating QCD 
coherence are below ~15%. Therefore, we evaluate Jo from Table [2] as 

e + e" annihilation : J e+e_ w (3.3 ± 0.5) GeV~ 2 . (16) 



8 The Gaussian parametrization of the distribution over the rapidity y seems to be a bit low at y > 1 (see 
Fig. 24 of [21]), but this affects the spectrum shown in Fig. 2 by less than ~ 5%. 

9 The data plotted in [5] refer to [tt + + 7r~], and we have to multiply the results by a factor 1/2. 
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The error in this result should include the uncertainty caused by the extrapolation ( 115]) over 
the energy E, which also allows to fit the data at very low momenta p ~ 0.05 GeV. 

Table 2: Soft limit Jo in the e + e~ collisions of the fits ({15"]) (leading term Ai, B{) from TASSO [33J 
to the single pion energy spectra \{it + + 7r~)/2] and the estimates, found using data from 
ARGUS [3]Zj and TPC/2 7 [31] . 



Exp 




Pmin 


Ax 


B x 


Otot 


h 




[GeV] 


[GeV] 


[nb/GeV 2 ] 


[GeV- 1 ] 


[nb] 


[GeV- 2 )] 


ARGUS 


10 


p > 0.05 








3.0 ±0.5 


TASSO 


14 


p > 0.3 


23.9 ±3.9 


5.25 ±0.43 


1.77 


3.5 ±0.6 


TASSO 


22 


p > 0.3 


8.0 ± 1.2 


4.70 ±0.32 


0.72 


3.4 ±0.5 


TASSO 


34 


p > 0.3 


3.7 ±0.6 


4.97 ±0.45 


0.30 


3.4 ±0.5 


TPC/27 


29 


p T ~ 0.05 








3.0 ±0.3 



It is also interesting to compare with DIS results in the current fragmentation region in the 
Breit frame, which should correspond to the quark fragmentation similar to e + e~ annihilation. 
The HI Collaboration at HERA [31] performed a study of this type at Q 2 varying from 12 to 
100 GeV 2 . Indeed, they found that the particle density at low momenta is nearly constant. 
The soft limit is obtained by extrapolation of the data in energy towards E = Qq = 270 MeV, 
which is used as an effective mass of charged particles. Extrapolating the data points at the 
lowest energy above E = 0.3 GeV to E = Q at < Q >= 19.6 GeV we arrive at J ~ 4 ± 0.5 
GeV" 2 , which is compatible with Eq. (fl6l) . It would be interesting to analyse the recent more 
detailed data from ZEUS [35] in a similar way. 

Finally, we can compare with the pp collision data, using the estimates (fill) and (fl2|) for 
the non-diffractive events. Then we obtain 

C/Io +e ~ » (1-8 ± 0.4) + (2.4 ± 0.5), (17) 

where the two numbers correspond to the thermal or exponential parametrization respectively. 
As discussed above, the evaluation based on the exponential fit may be considered as an upper 
limit, which would apply for a distribution without accounting for a mass effect. The uncer- 
tainty in (j!7p will be largely reduced by more precise data on the spectrum at very small pt in 
pp collisions. 

The result ffTTj) agrees well with the expectation from the different primary quark and gluon 
sources in (jSJ). 

The approximate energy independence of the quantity 1$ in both collision processes, as well 
as the obtained ratio of the Iq values in the pp and e + e~ reactions, which is about 2, are remark- 
able. This is a serious argument in favour of the relevance of the elementary bremsstrahlung 
process, also for the soft pp interactions. 
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4 Nucleus-Nucleus interaction 



4.1 Spectra at low transverse momenta 

In the case of AA scattering one would naively expect, that the soft particle density is equal 
to that in the pp-collisions times the mean number of nucleon-nucleon collisions, N co u. Such 
an estimate should be valid for the point-like interactions, for example, for the events with the 
large particle production, if we neglect the energy losses in the nuclear medium. However, 
at low momentum transfer, p?, the coherence effects should reduce the particle production, 
in particular, due to the destructive interference between the different amplitudes within the 
space domain of the size of ~ 1/pt • Thus, we can expect the limiting behaviour when the soft 
particle wave-length, 1/pr, becomes comparable with the coherence range, at most with the 
nuclear radius ta (for the central AA-collisions) or with the size of the region, where two nuclei 
overlap, that is, for instance, for the AuAu collisions pt < ^/ta ~ 30 MeV. Inspection of the 
RHIC data [SB] - [39] shows, that, indeed, the ratio 

r n cM 1 dN AB /dp T , , 

N coll dN pp /dp T ' 

decreases with decreasing transverse momentum of the secondaries. The ratio Rab u fllBD 
allows to compare the particle density for the nucleus-nucleus collisions with that in the proton- 
proton interactions times the mean number of nucleon-nucleon collisions N co u, which is calcu- 
lated using the Glauber model. Another way used to present the RHIC data is to replace the 
number of collisions N co u by the number of nucleons participating in the interaction, N part . 
Since in the pp-collisions just two protons scatter, the normalization factor is N part /2 rather 
than Ncou. 

Similar to the limiting distributions, corresponding to the soft bremsstrahlung in the e + e~ 
and pp processes, we, first of all, expect the energy-independent behaviour of the soft particle 
densities. The existing data confirm, that at low transverse momenta, pt < 0.3 GeV, the energy 
dependence of the ratios Rab is, indeed, quite flat. At so low p^s the ratios Rab, measured by 
PHOBOS [39] at y/s = 62.4 and 200 GeV, largely coincide for all studied centrality ranges (see 
Fig. 8 and Fig. 32 in [39] for the normalizations with N co u and N part respectively). At lower 
energies y/s^N = 5 and 17 GeV the results were obtained by PHENIX [37] (see Fig. 47), and 
a comparison is possible for pt > 0.4 GeV. A convergence of the pt spectra towards low px at 
the two energies can be seen, although not yet a coincidence in the observed pr-range. 

Let us now turn to the limiting density of the soft particle production. The data at the 
lowest p T « 0.03 GeV were collected by PHOBOS [39] (Fig. 3), the data with p T > 0.2 GeV 
were obtained by PHENIX [37] (Fig. 11) and by STAR [27] (Fig. 46, Tab. 33) in the Au+Au 
collisions at 200 GeV. If we extrapolate the invariant pion spectra from PHENIX and STAR 
by the exponential form at pr < 0.5 GeV we arrive at l AuAu « 1270 GeV~ 2 to be compared 
with Iq P « 7.8 GeV -2 , found with the same procedure for the pp collisions [27J. However, 
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the PHOBOS data at very low p? indicate a flattening of the distribution according to the 
functional form © with T = 0.229 GeV. 

Fitting the PHOBOS and STAR AuAu data for pions at < 0.5 GeV using this parametriza- 
tion, we obtain for the mean pion density at pt = in the average 

central AuAu collisions (pions) : I ~ (950 ± 100) GeV" 2 , (19) 

which is about 30% lower than the result of exponential parametrization. Assuming a similar fit 
for the low pt pp STAR data, we arrive at Iq P ~ 5.9 GeV -2 with T = 0.182 GeV, therefore, the 
limiting density for the nuclear collisions is more than 100 times larger than the corresponding 
value in the pp interactions. On the other hand, in this case the Glauber model calculation [39J 
gives N co u = 1040 and N part /2 = 172 (±15%). Clearly, the soft particle density is about 10 
times lower than that expected in the case of independent (incoherent) collisions. Indeed, from 
the extrapolation of the pt spectra in the AuAu and pp collisions by either an exponential or 
the thermal function (Q, as presented above, we find the ratio 

Io A /I™ « 160 ± 17, (20) 
which agrees with the calculated N part /2, and, therefore, 

p T ^0: R%? rt ^l and l£ uAu « ^ J PP . (21) 

Here R^ rt is defined as in (TlBl) . but with N part /2 as a normalizing quantity. 

A detail experimental study was performed by PHOBOS [39] (Fig. 32) for different ranges 
of centrality, that is, the number of participants, at the two energies 62.4 GeV and 200 GeV. 
Remarkably, the ratio Rj^ Tt approaches about unity in the soft limit for all selections of 
centralities or the N part parameter. STAR [41] measured the ratio R c Xa with high precision for 
Pt > 0.5 GeV. It is found that this quantity decreases below px < 2 GeV and the extrapolation 
below the measured values of px suggests the limiting behaviour f[2"Tj) as well. 

That is, the data are consistent with the model where each pair of participating nucleons 
produces its own number of secondaries independently of the number of collisions each nucleon 
is participated in. This is related to the idea, that the overall characteristics of bulk particle 
production depend only on the number of 'wounded nucleons' [40] and not on the number of 
rescatterings. 

4.2 Antenna pattern in nuclear collisions 

Within the bremsstrahlung scenario this result suggests that there is a coherent particle pro- 
duction over the range of nucleon size, but different nucleons are separated in space-time so that 
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their contributions remain incoherent This is in agreement with the observation that no spe- 
cific structure reveals itself in the nuclear collisions at very small px ~ 30 MeV corresponding 
to the nuclear size (see, for example, Fig. 3 in [35] with pions in this range). 

In spite of the multiple re-interactions of the same nucleon in the nucleus which proceed 
through the corresponding number of gluon exchanges in the minimal model for the pp collisions 
considered above, we still expect the limiting behaviour. The low pt particles emitted in this 
process interact coherently with all t-channel gluons and, actually, probe the overall colour 
flow originated by this (t-channel) system of gluons. On the other hand, as we discussed in 
Sect. 3.1, even a large number of gluons produces dominantly the same octet colour flow as 
in the pp inelastic collision (first, the octet exchange gives the dominant Leading Logarithmic 
contribution, both in the DGLAP and in the BFKL kinematics, next - interacting with the 
individual quark or with the colourless qq- or 'quark-diquark'-system we can transfer only the 
colour octet quantum number). Examples of diagrams contributing in this minimal model to 
the production of the low pt particles in the pp and pA collisions are shown in Fig. [3J In the 
AA collisions re-interactions of nucleons may occur as in the pA collisions. In these depicted 
processes the colour flow and the associated soft particle production are determined by the 
number of participating nucleons irrespectively of the number of rescatterings, in agreement 
with the phenomenological result f[2"Tj) . 

4.3 Universality of particle ratios at low px 

Having in mind this universal mechanism for soft particle production we may ask whether there 
are any consequences for the relative rates of different hadron yields in the same limit. If the 
particles at low momenta in the hadronic and nuclear collisions are related to the universal 
bremsstrahlung from the (incoherent) superposition of primary colour octet charges, then, in 
this limit, the particle composition should be the same in the different processes. To the extent, 
that the dominant underlying mechanism in the pp collisions is the gluon exchange between 
the quark constituents of the proton with the initial and final bremsstrahlung, the e + e~ and 
hadronic data on the particle ratios should approach each other as well. 

A similarity of particle ratios K/n and p/ir in the e + e~ and pp reactions at pr < 0.5 GeV has 
been indeed noted already some time ago by the TPC collaboration [3TJ- In this measurement 
the transverse momentum pt for the particles in e + e~ collisions was defined with respect to 
the sphericity axis. 

The pt dependence of the particle ratios for several hadronic collisions have been compared 
by PHENIX [12]. While at the large pr > 2 GeV the ratios p/ir and K/ir tend to approach large 
values ~ 1 in the central AuAu collisions, these ratios are reduced for non-central and minimum 
bias pp collisions. Remarkably, these ratios converge for the different processes towards lower 
Pt < 1 GeV. In Fig. H] we collect data in the low px region on the ratio K~/tt~ from the 

10 The coherence is destroyed by the presence of different recoil nucleons. 
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Figure 3: Diagrams, contributing to the pp and pA collisions in the minimal model for soft 
particle production: (a) In the pp collisions the exchanged gluon interacts with the colour 
triplet constituents q or qq in the proton to form an outgoing colour octet system; (b,c) In the 
pA collisions the proton can rescatter inside the nucleus and then forms a colour octet system 
again. This implies that the multiple gluon exchange acts as a single gluon exchange in the 
particle production. In this example: N co u = 3, N part /2 = 2. 

e + e~, pp and AA interactions. As one can see, the particle ratios, indeed, approach each other 
towards low px < 0.4 GeV pointing towards a dominance of multiple qq dipole radiation in all 
processes. 



5 Conclusions 



An issue of universality of particle production in various collision processes predated QCD. The 
perturbative dynamics based on the QCD gluon bremsstrahlung suggests a particular type of 
universality, which, however, allows only for the predictions for the soft limit px — > 0. This is 
the case, when the lowest order Born term with its energy independence and elementary colour 
factors dominates, whereas for the larger p^s the (non-universal) higher order contributions 
take over. The predictions from this approach are confirmed with high accuracy in the e + e~ 
measurements, where the energy independence of the soft particle yield is observed, and its 
intensity in quark and gluon jets is in agreement with the expectation from corresponding 
colour factors Cp and Ca- 
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Figure 4: Convergence of particle ratios K /n towards small Pt for various processes: e + e 
annihilation (TPC (31] data, px with respect to sphericity axis), pp (minimum bias) and central 
(0-5%) AuAu collisions (STAR [27] and PHENIX [32] Collaborations). 

In this paper we discuss the production of low momentum particles in hadronic collisions. 
Surprisingly, the soft particle production density in the pp collisions is practically constant over 
the large energy range up to 1800 GeV. Moreover, the approximate energy independence holds 
also in the nuclear collisions at RHIC energies. The colour factors related to the primary in- 
teraction process reveal themselves also in the comparison between the e + e~ and pp collisions, 
according to a minimal model with gluon exchange in the case of hadronic processes. Such a 
comparison would strongly benefit from more precise measurements of the low pt particle spec- 
tra in pp collisions, which could reveal the expected deviations from the exponential behaviour. 
In the nuclear collisions the perturbative expectations, based on the coherent bremsstrahlung 
in a multi-gluon exchange process lead to the energy independence of the soft particle yield. 
Moreover, the magnitude should scale with the number of participant nucleons, N part . 

As a result, the yield of the low px particles can be directly related to the soft QCD 
gluon bremsstrahlung from the colour charges, which are created in the primary hard or semi- 
hard interaction, and is proportional to Cf, Ca and Np ° rt Ca for e + e~, pp and AA collisions 
respectively. This concept of universality is also supported by the observed convergence of the 
7r : K : p particle ratios in all these processes in the pr — > limit. 

From the first sight, this seems to contradict the idea of thermalisation in the particle 
production. However, we recall, that according to the space-time picture of particle production, 
the central soft hadrons in the nucleon-nucleon c.m.s. system are those, which are formed first, 
and these soft particles with their large wave-length probe only the colour charges participating 
in the primary local interaction at a transverse size of ~ 1 f. The fast particles are formed 
later on, and leave the interaction region without re-scattering on the slow hadrons which stay 
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in the interaction region. Therefore, the slow particles may be not in a thermo-equilibrium 
with the whole hadronic system. At the very early stages following the primary interaction 
no re-interactions of produced hadrons from the different nucleons will take place, while the 
re-interactions of one nucleon with other nucleons in the nucleus should not be accounted for 
because of the coherent nature of soft particle emission. Then, the universality of this soft 
particle production from the primary colour sources becomes a plausible phenomenon, since 
it basically reflects the universality of soft particle production from the individual isolated 
quark-antiquark dipoles. 

It will be very interesting to extend such measurements on the limiting soft particle pro- 
duction in pp and AA collisions to the higher energies at the LHC and to see whether any new 
incoherent sources appear. Such measurements could set a critical benchmark for the models 
of multiparticle production in hadronic interactions. 
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